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Abstract 

We obtain a rigorous a priori upper and lower bounds to the exact period of the 
celebrated Rayleigh stretched-string differential equation. 

1 Introduction 

In his monumental treatise The Theory of Sound (see pages 45-46 [2J) Lord Rayleigh 
proposes the following problem: 

A horizontal elastic steel wire having length 2Lq and spring constant a is stretched 
to length 2L. A body, having a mass m much greater than the mass of the wire, 
is tied to the center of the wire and is put in motion on a vertical line through 
the center of the wire. The only forces on the body are produced by tension in the 
wire (there is no gravity, no damping). It is required to study the motion of the 
body. 

(We have used Agnew's formulation JT\, p. 17-18.) Let y{t) = yhe the vertical displacement 
of the mass m in time t, where y is positive, 0, or negative according as the body is above, 
or at, or below the wire in equilibrium or neutral position. 

By the law of HoOKE, the magnitude of the elastic force, or the tension, T{y), exerted 
by each half of the stretched string is equal to: 

^, , stretch 
T{y) = a- 



a ■ 



original length 



and therefore the vertical component of the total force, exerted by both halves, which is the 
only component acting to move the mass m is given by 

vertical component of T(y) = —2a- ^ ~ ~ — ^ ■> (1-1) 

Lq a/ L2 + y2 



where we have employed our sign conventions. 

But, by Newton's second law of motion, that same force, with opposite sign, is given 
bym||. 

Therefore, equating the two forces, and bringing the expression for the elastic force on 
the same side as the force given by Newton's second law, we find that the algebraic sum 
of the forces is zero and dividing by m, we conclude that the equation of motion of the mass 
m is given by: 



-TT H y = 

dt^ m Lo^/L'^ + y^ 
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This differential equation fll.2p is nonlinear and very complicated, and we have never seen 
an exact treatment of it in the literature. 

Rayleigh, himself, simplifies the equation by the following physical reasoning; 

"The tension of the string in the position of equilibrium depends on the 
amount of the stretching to which it has been subjected. In any other position 
the tension is greater; but we limit ourselves to the case of vibrations so small that 
the additional stretching is a negligible fraction of the whole. On this condition, 
the tension may be treated as a constant." 

(Our italics.) We just saw that the magnitude of the tension, T{y), (or the force) from 
each half of the string on m is given by 



T{y) = a 
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provided the stretching is well within elastic limits. So, Rayleigh effectively takes y = in 
(11. 3p so that the equation (II. 2p becomes 



(Py 2a L — Lq 
dt"^ m LqL 



y = o, 
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or, using the notation (11.31) (with y = 0, i.e., T := T(0) = cr— ^), 

2T 

y + —y = 
mL 



;i.5) 



where the "dots" refer to time derivatives. If we assume that the initial displacement at time 
t = is yo, then the solution to (11.51) is 
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which is simple harmonic motion with period P given by: 

271 



P 
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Rayleigh does not discuss how close his approximate period fll.7p is to the exact period 
P. Indeed, we have been unable to to find any error analysis of Rayleigh's solution in the 
literature. Yet, such a famous and classical differential equation merits an investigation into 
the accuracy of its approximate solution. 

Therefore we offer the following theorem to fill this gap. 

Theorem 1. Let P be the period of oscilation of the mass m. Then Rayleigh 's approximate 
period (11. 7p overestimates the true period. Indeed, the following inequalities are valid: 




Our lower bound shows the dependence of the period on the initial displacement, Uq, and 
tends to the upper bound as Uq tends to 0. 



2 The Formula for the Period 

The differential equation (11.21) does not have the variable t nor the first derivative, y , ap- 
pearing explicity. Therefore, we can transform it into a separable equation for the velocity, 
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(2.1) 
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dv dv dy 
dt dy dt dy 

Substituting the right-hand side of f l2.2p into ( 11. 2p . transposing the expression for the tension 
to the right-hand side, separating the variables and forming the indefinite integral of both 
sides we obtain 



(2.3) 



where C is the constant of integration. At time t = 0, the conditions are y{0) = yo and 
f (0) = 0. Therefore, evaluating the constant C, recalling that v = taking the square 
root of both sides, and using a negative sign because the displacement, y{t), is a decreasing 
function during the first quarter oscilation, we obtain: 



dy 
dt 



\ 



m 




(2.4) 



Multiplying the quarter-period by 4, we obtain the following formula for the full period, P, 
of oscilation of the mass: 
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Theorem 2. The true period, P, of oscillation of the mass m is given by: 



P = 4 




yo 



dy. 



This integral is very comphcated for exact computation. 
Indeed, if we make the change of variable 



z^:=L^ + y^ zl:=L' + yl 



we obtain the formula 



P = 4W- 

2a ,L 



dz. 
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If we write 



as 
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and 



then 



Pf^) :=4W — 
2a 



dz 



Lj^,-a){z-b){z-c){z-d) 



P = P{zo)-P{L). 



(2.5) 



(2.6) 



(2.7) 
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The integral in 02.91) is a very complicated algebraic combination of an elliptic integral of 
the third kind and an elliptic integral of the first kind, and we will not pursue this line of 
inquiry any further. 

Instead, we turn our attention to the upper and lower bounds for the true period, P. 



3 Upper and Lower Bounds for the True Period 

We begin with the exact formula (13.11) for the true period, P: 




(3.1) 



Taking y = yo = in the sum a/XM-"^ + ^/iX^-y^ only makes the integrand larger 
and we conclude that 
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P <4 




yo 



dy 



T 2 



2tt 



and ti;e /iave proved the validity of the upper bound: 




(3.2) 



For the lower bound, we proceed similarly, and we start by observing that if we take 
y = yo in the sum L? ^ y"^ + ^J + y"^ the integrand becomes smaller and we conclude 
that 
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where, twice, we have used the inequality 



L2 + y2 ^ L + 



2L' 



Therefore, we have proved the validity of the lower hound: 




(3.3) 



This completes the proof of the theorem 11.81 
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Corollary 1. If P denotes Rayleigh's approximate period {\1.7\i . then, the relative error 

P-P 



R :-- 



P 



in the approximation P ^ P, satisfies the inequality 



vim 
"4TLn 



^ i? ^ 



(3.4) 



This shows that the relative error is proportional to the square of the initial displacement 
and tends to zero with a "quadratic" rate of convergence. 
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